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Relaxation of a Single Knotted Ring Polymer 

Shinya Saka* and Hiroshi Takano 
Faculty of Science and Technology, Keio University, Yokohama 223-8522, Japan 

The relaxation of a single knotted ring polymer is studied by Brownian dynamics sim- 

O ' . . 

' ulations. The relaxation rate X q for the wave number q is estimated by the least square 

fit of the equilibrium time-displaced correlation function C q (t) = TV -1 J2i Ylj ' 

Rj(0)) exp[i2Tr q(j — i)/N] to a double exponential decay at long times. Here, N is the number 

of segments of a ring polymer and Ri denotes the position of the zth segment relative to 

the center of mass of the polymer. The relaxation rate distribution of a single ring polymer 

with the trivial knot appears to behave as X q ~ A(l/N) x for q = 1 and X q ~ A'(q/N) x for 

q > 1, where x ~ 2.10, x' ~ 2.17 and A < A'. These exponents are similar to that found 

for a linear polymer chain. The topological effect appears as the separation of the power law 

dependences for q = 1 and q > 1, which does not appear for a linear polymer chain. In the 

case of the trefoil knot, the relaxation rate distribution appears to behave as X q ~ A(l/N) x 

for q = 1 and X q ~ A'(q/N) x ' for q = 2 and 3, where x ~ 2.61, x' ~ 2.02 and A > A'. 

The wave number q of the slowest relaxation rate X q for each TV is given by q = 2 for small 

values of N, while it is given by q — 1 for large values of AT. This crossover corresponds to 

the change of the structure of the ring polymer caused by the localization of the knotted 

^ • part to a part of the ring polymer. 

KEYWORDS: knot, ring polymer, single polymer, relaxation modes, relaxation rates, Brownian 
dynamics simulations, topological effects 
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1. Introduction 



X 

5-i 



The effects of topological constraints caused by the entanglement of polymers on the 
properties of polymer systems have attracted much attention. 1-21 One of the famous theories 



on the topological effects is the reptation theory 1 ' 2 of the dynamics of concentrated polymer 
systems, where polymers are entangled each other. In this theory, the topological constraints 
on a polymer chain caused by the surrounding polymer chains are replaced by an effective 
tube confining the polymer chain. The predictions of the theory with corrections due to the 
finiteness of polymer length 2 agree with the experimental results and have been confirmed 
by simulations. 3 ' 4 In concentrated polymer systems, the entanglement of polymers changes 
dynamically and the topological constraints change accordingly. In contrast, in the case of 
a single knotted ring polymer, which is one of the self-entangled systems, the topological 
constraints are determined by the knot type and do not change with time. Therefore, a single 
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ring polymer system can be considered as an ideal system for the study of topological effects. 
The investigation of the effects of knots in this system is expected to provide a basis for further 
understanding of the topological effects in polymer systems. 

The interest in the properties of knotted ring polymers has been growing in recent 
years. 5-21 Especially, the topological effects on the static properties of ring polymers have 
been well studied. 6 " 15 It was conjectured by des Cloizeaux that topological constraints make 
ring polymers with long chains swell. 6 Although this swelling effect caused by effective repul- 
sion among segments is supported by a computer simulation 7 and a scaling analysis, 8 there 
have been studies 9-13 predicting that the effect vanishes in the long chain limit. In contrast 
to the static properties, the topological effects on the dynamics are less well studied. 16-18 The 
equilibrium relaxation of single knotted ring polymers has been studied by measuring the 
time autocorrelation function of the radius of gyration of a knotted ring polymer and the long 
relaxation time that does not appear for the trivial knot has been found. 16 The relaxation 
time has the dependence on the number of the essential crossings C. 17 The nonequilibrium 
relaxation after cutting one bond of a knotted ring polymer has been studied and the dis- 
tribution of the nonequilibrium relaxation time of the radius of gyration has been found to 
show different behavior for different knot groups. 5,18 Experimentally, knotted circular DNA 
molecules have been studied by using gel electrophoresis, electron microscopy, and so on. 19-21 

In the case of linear polymers, the relaxation phenomena have been studied systematically 
in terms of the relaxation modes and rates. 3 ' 22-2 The relaxation modes and rates are given as 
left eigenfunctions and eigenvalues of the time-evolution operator of the master equation of the 
system, respectively. 22-24 The equilibrium time correlation functions of the relaxation modes 
satisfy (X p (t)X q (0)) oc 5 Ptq exp(— X p t), where X p and X p denote the pth relaxation mode and its 
relaxation rate, respectively. For single linear polymers represented by the Rouse model, which 
has no excluded volume interaction and no hydrodynamic interaction, the relaxation modes 
are identical to the Rouse modes, which have been playing an important role in the theory of 
the polymer dynamics. 2 ' 28 ' 29 For single linear polymers with the excluded volume interaction 
and without the hydrodynamic interaction, the relaxation modes and rates are estimated by 
solving a generalized eigenvalue problem Ylf=i Cj,j(io + T )fp,j = exp(— A p r) Ylf=i Ci,j{to)fp,j 
under the orthonormal condition Y^iLi YljLi fp,iCi,j{to)fq,j = fip,q, 23 ' 24 where Cjj(i) denotes 
the equilibrium time-displaced correlation function ^(Ri(t) ■ Rj(0)} of the position of the 
ith segment relative to the center of mass of the polymer R4 and that of the jth segment 
and N denotes the number of the segments of the polymer. The pth relaxation mode and 
the corresponding relaxation rate are given by f p ^ and X p , respectively. The behavior of the 
contribution gi :P of the pth slowest relaxation mode to Ri is similar to that of the Rouse 
mode, gi tP oc cos [irp(i — |)/AT], and the corresponding relaxation rate X p behaves as X p oc 
(p/A) 2l/+1 . 2 ' 23 ' 24 Here, v ~ 0.588 is the exponent for the power law dependence of the size of 
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a single linear polymer with the excluded volume interaction on the number of the segments 
of the polymer. 

The purpose of the present paper is to examine the effects of the topological constraints 
on the relaxation of single ring polymers by carrying out an analysis similar to that which has 
been done for linear single polymers. Brownian dynamics simulations of single ring polymers 
with the trivial knot and the trefoil knot are performed, where the excluded volume interaction 
is taken into account and the hydrodynamic interaction is neglected, and the distribution of 
relaxation rates is estimated. 

In §2, a model used in the present paper and the method for the estimation of the relaxation 
rates are explained. The results of the simulations are presented in §3. Summary and discussion 
are given in the last section. 

2. Model and Relaxation Rates 

In order to study a single knotted ring polymer in good solvent, Brownian dynamics 
simulations of a bead-spring model are performed. The dynamics of the ith segment of a 
single ring polymer with N segments is described by the Langevin equation 

dn(t) __ldV({ rj }) 

c 



+ Wi (t). 



(1) 



dt C dri 

Here, ri{t) denotes the position of the ith segment at time t and £ is the friction constant. 
The random force (Wi(t) acting on the ith segment is a Gaussian white stochastic process 
satisfying (wi t0l (t)) = and 



(2) 



where Wi j(X , fee and T denote the a-component of Wi, the Boltzmann constant and the tem- 
perature of the system, respectively. The potential V({rj}) = F(n,--- ,rjv) describes the 
interaction between the segments. In eq. (1), the hydrodynamics interaction is neglected. 
In the present paper, we use the potential given by 24, 30-32 

V ({rj}) = £E V R (\n - r s \) + V A (\r i+1 - ri \) (3) 
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Here, Vr is the repulsive part of the Lennard-Jones potential and represents the excluded 
volume interaction between all the segments. The potential V\ is called a finitely extensible 
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nonlinear elastic (FENE) potential and represents the attractive interaction between neigh- 
boring segments along the ring polymer. 32 In the last summation of the right hand side of 
eq. (3), rjy + i = r\ because the Nth. segment is connected to the first segment along the ring 
polymer. 

As mentioned in the previous section, the relaxation modes and rates are given as left 
eigenfunctions and eigenvalues of the time-evolution operator of the master equation of 
the system, respectively. 22-24 In the previous studies on the relaxation modes and rates 
of single linear polymers, a trial function for the pth relaxation mode is chosen to be 
X p (Q) = ^2iLi fp,iRi{to/2;Q), where Q represents a state of the system and Ri(to/2;Q) 
denotes the expectation value of Ri after a period to/2 starting from a state Q. The quantity 
Ri denotes the position of the ith segment relative to the center of mass of the polymer: 
Ri = f i — t c with r c = X^^i r i- For this trial function, a variational problem equivalent to 
the eigenvalue problem for the time-evolution operator leads to a generalized eigenvalue prob- 
lem ~}2!j = iCi y j(tQ +r)f p j = exp(— A p r) Ylj=i Ci,j(to)fp,j mentioned in the previous section, 
where Cij(t) = |(J2j(t) ■ Rj(0)) is the equilibrium time-displaced correlation function. This 
analysis is considered to extract the slow relaxation modes contained in the quantities {Ri} 
and to give better results as to becomes larger, because the contribution of faster relaxation 
modes contained in Ri(to/2;Q) becomes smaller. 

In the case of single ring polymers, the eigenfunctions f p j of the above-mentioned gener- 
alized eigenvalue problem are known because of the translational invariance of the segment 
number along the ring polymer: Cij(t) = Ci + ij + i(t), where the subscripts representing the seg- 
ment numbers are considered modulo N. The eigenfunctions are given by f q j oc exp (i2irqj/N) 
with q = l,--- , N - l. 33 Note that the wave number q = is not included for it gives a mean- 
ingless relaxation mode ~}2^ = iRj{to/2;Q) = with an eigenvalue exp(— A g r) = 0, because 
J2f=i Rj = an d Ylj=iCi,j{t) = 0. For the eigenfunction with the wave number q, the 
generalized eigenvalue problem is reduced to 

C q (to + r) = exp(-X q T)C q (t ), (6) 

where C q (t) denotes the Fourier transform of the correlation matrix given by 

N-i / 2 v 

C q (t) = £C M+J (t)exp (i—ql) (7) 

= £ C M+l(*)co6 l—ql • (8) 
1=0 ^ ' 

Note that Cj does not depend on i and is identical to Ci t i-i(t). The equilibrium time- 

displaced correlation function of the Fourier transform R q of Ri, which is defined by 

' (9) 



1 - /2tt 
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is also given by C q {t) as 



C q (t) = -(R q (t)-R- q (0)). (10) 



1 

3' 

In the following, we only consider C q (t) with q = 1,2, •■ ■ , [N/2\, since C q (t) = Cjv_ g (i) as 
can be seen from eq. (8). Here, [x\ denotes the floor function of a real number x. Equation (6) 
is equivalent to the estimation of the relaxation rate X q for the wave number q from the values 
of C q (t) at times to and to + r by assuming a single exponential decay of C q (t) at these times. 
In fact, C q (t), which is the autocorrelation function of the Fourier transform of R4, can be 

expressed as C q (t) = exp(— A^t) + a q 2 ^ exp(— A^t) H , where < A^ < X q 2 ^ < ■ ■ ■ and 

a q k ^ > for k = 1, 2, ■ • • , 22 ' 23 Thus, the estimation of the slowest relaxation rate A^ by eq. 
(6) becomes better as to becomes larger, as mentioned before. In this paper, we estimate A^ , 
the slowest relaxation rate in C q (t) for the wave number q, by fitting several tens of values of 
C q (t) at relatively long times estimated from the simulations to a double exponential decay 
form 

C q (t) = aW exp(-A«t) + a< 2 > exp(-A( 2 )t), (11) 

because this method is more tolerant of statistical errors in the estimated values of C q {t) than 
using eq. (6). 

In the linearization approximation of Doi and Edwards, 2 ' 23 the relaxation rate \ p of the 
pth Rouse mode X p = (2/N) 1/2 Ri cos [pir(i - \)/N] of a single linear polymer with 
the excluded volume interaction is estimated from the static correlation (X p ■ X p ) as X p LA ~^ = 
(fceT/C) {\{X P ■ X p )) \ In the case of single ring polymers, the relaxation rate X q for the 
wave number q is given by the linearization approximation as 

Af A > = ^QA-fl- ff >)" 1 = ^(0)- 1 . (12) 

These formulae are derived as follows. 2 In the linearization approximation, the equation of 
motion of each mode x is approximated by the linear Langevin equation x = —\x + w, where 
the Gaussian white random force w(t) satisfies (w(t)) = and (w(t)w(t')) = 2D5(t — t'). On 
the basis of the relation (x 2 ) = D/X, which holds for this Langevin equation, the relaxation 
rate A is estimated from the mean square amplitude of the mode (x 2 ) as A = D/(x 2 }. Thus, 
the relaxation rate given by the linearization approximation becomes small as the mean square 
amplitude of the mode vector, X p or R q , becomes large. The mode vector X p for the linear 
polymer roughly corresponds to the end-to-end vector of a partial chain with a length of 
N/p segments and R q for the ring polymer roughly corresponds to the end-to-end vector of 
a partial chain with a length of N/(2q) segments. For a single linear polymer, the results 
of the linearization approximation are found to be consistent with those of the Monte Carlo 
simulations 23 and the molecular dynamics simulations. 24 In this paper, the relaxation rates 
Aq LA ^ given by the linearization approximation for a single ring polymer are estimated by 
using eq. (12) and compared with the relaxation rates Xq 1 ^. 
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3. Results of Simulations 

Brownian dynamics simulations of the model described in the previous section are 
performed for N = 24, 36, 48, 72, 96, 144 and 192. The following parameters 24 ' 31 are 
used: k-^T/e = 1, ka 2 /e = 30 and Rq/ct = 1.5. The Euler algorithm with a time step 
At = 1.4 x 10 -4 7(7 2 /e is employed for a numerical integration of the equation of motion 

(I) . Hereafter, we set a = 1, 7 = 1 and e = 1. 

The correlation function C q (t) is calculated from Ci : j(t) = ^(Ri(t) ■ Rj(0)) estimated by 
the simulations by using Eq. (8). The equilibrium average (Ri(t) ■ Rj(0)) is estimated as the 
time average over M c configurations after the initial Mi configurations, which are discarded 
for the equilibration: 

Mi+M c 

(Ri{t) ■ RjiO)) = w Yl Ri(nAT).Rj{nAT-t), (13) 

c n=Afi+l 

where the configurations are taken from a simulation at intervals of time AT. In the present 
simulations, M u M c and AT are chosen as M^AT = 10r, M C AT = lOOOr and AT = 2 x 1(T 3 t, 
where r is about the longest relaxation time l/A^i of the ring polymer with the trivial knot 
for each N. 

The relaxation rate Aq 1 '* are estimated by the weighted least square fit of C q (t) to eq. 

(II) . The range of t used for each fit satisfies the conditions 0.25 < C q (t) / C q (0) < 0.95 and 
5Cq(t)/C q (t) < 0.1, where SC q (t) denotes the statistical error of the estimate of C q (t). 

Figures 1(a) and 1(b) show log-log plots of X q ^ versus q/N and X q LA ^ versus q/N, respec- 
tively, for the ring polymer with the trivial knot, where the number of segments N = 24, 
36, 48, 72, 96, 144 and 192. The behaviors of X q LA ^ in Fig. 1(b) are qualitatively similar to 
those of Xg 1 ^ in Fig. 1(a), although each value of X q LA ^ does not quantitatively agree with the 
corresponding value of Xq 1 ^, in contrast to the case of a single linear polymer. 23,24 In each of 
Figs. 1(a) and 1(b), the data points for q = 1 and those for q > 1 seem to fall on two different 
straight lines at small values of q/N. This suggests the power law behaviors 

(*) 



.(*) / 1 



A iIn) fOT 9 = 1 . 



„(*) 



(14) 

41 (ir fa o 1 - 

where the superscript (*) denotes (1) or (LA). The amplitudes and exponents in eq. (14) are 

(*) 

estimated by the least square fit of the data points to a straight line in the log-log plot: A\ ' 
and x^ 1 are estimated from the data points for q = 1 with iV = 48-192; A^*\ and x>\ are 
estimated from the data points for q > 1 with q/N < 0.05. The estimated parameters for the 
relaxation rates A^ 1} are given by ~ 24.0, ~ 2.10, ~ 50.2 and ~ 2.17. The 
exponents x^p and x>\ are similar to that for a linear polymer chain 2u + 1 ~ 2.18. 2,23 ' 24 The 
separation of the power law dependences for q = 1 and q > 1 is attributed to the topological 
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the trivial knot, where N = 24, 36, 48, 72, 96, 144 and 192. For each value of N, the solid symbols 



represent the relaxation rates for q = 1 and the open symbols represent those for q > 1. In each 
figure, the solid line represents the result of the least square fit of the data points for q = 1 with 
N = 48-192 to a straight line and the broken line represents that for q > 1 with q/N < 0.05. 



constraints, because no such separation appears for a single linear polymer 23 ' 24 and an ideal 
single ring polymer, 33 which have no excluded volume interaction. A similar separation of the 
power law dependences of the relaxation rates has been observed for a single linear polymer 
trapped in an array of obstacles in two dimensions. 27 The parameters for Aq LA ^ are estimated 

as Af A) ~ 60.9, xf A) - 2 - 19 > A >\ ] - 214 and x >^ - 2 - 39 - Note that the relations 

(*) (*) (*) (*) (1) (LA) 

x\ < x>{ and A\ < hold for both \\ and \ q .It follows from the relations that the 

two straight lines in the log-log plot, which correspond to the power law behaviors for q = 1 
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0.01 0.1 
q/N 



Fig. 2. Log-log plots of (a) A, ' versus q/N and (b) \ q ' versus q/N for the single ring polymer 
with the trefoil knot, where N = 24, 36, 48, 72, 96, 144 and 192. For each value of N, the solid 
symbols represent the relaxation rates for q = 1 and the open symbols represent those for q > 1. 
In each figure, the solid line represents the result of the least square fit of data points for q = 1. 
The broken line represents that for q = 2 and 3 in (a) and that for q = 2 in (b). The dash-dotted 
line in (b) represents that for q = 3. The data points for N = 48-192 are used for each fit. 



and q > 1, intersects at a small value of q/N. This suggests that the two power law behaviors 
for q = 1 and q > 1 are only apparent and eventually merge into a single power law behavior 
Aq*^ oc (q/N) x with X = 2u + 1 at sufficiently small values of q/N. 

In Figures 2(a) and 2(b), the g/iV-dependences of A^ and A^ LA ' ) are shown, respectively, 
for the ring polymer with the trefoil knot by a log-log plot, where the number of segments 
N = 24, 36, 48, 72, 96, 144 and 192. In the same way as in Fig. 1, although values of Af A) 
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in Fig. 2(b) do not quantitatively agree with those of in Fig. 2(a), the behaviors of Ag LA ^ 
are qualitatively similar to those of In Fig. 2(a), the data points for q = 1 and those for 
q = 2 and 3 seem to fall on two different straight lines for small values of q/N, while the data 
points for q = 1, those for q = 2 and those for q = 3 seem to fall on three different straight 
lines in Fig. 2(b). Therefore, we assume the power law behaviors 

A«^ 4*) (|) 4 *' for g = l, 2 and 3 (15) 

(*) (*) 

and estimate A\ and x q by the least square fit of the data points for each value of q with 
N = 48-192 to a straight line in the log-log plot. The parameters for A^ are estimated as 
A^ ~ 591, 4 1} - 2 -61, A^ ] ~ 7.54, ~ 2.02, A^ ~ 7.73 and x^ ~ 2.04. As expected, 
the relations ~ and 3 hold, which suggests that the relaxation rates Aq 

for q = 2 and 3 follow the same power law 

A«^ 4;](|) 4!1 for g = 2 and 3. (16) 

The least square fit of the data points for q = 2 and 3 with N = 48-192 gives A^\ c± 7.38 
and 45 - 2.02. The parameters for xf A) are estimated as Af A) ~ 2.97 x 10 3 , xf K) ~ 2.77, 
4 LA) ~ 16.8, xf K) ~ 1.95, Af A) ~ 58.8 and x^ LA) ~ 2.18. The estimated values of A 2 LA) , 
x^ A \ A^ A ^ and x 3 LA ^ lead to the intersection of the two straight lines for q = 2 and 3 at 
(q/N) ~ 4.3 x 10~ 3 , which suggests that the two apparent power law behaviors merge into a 
single power law behavior at small values of q/N in the same way as A^ with q = 2 and 3 
follows the power law (16). In Fig. 2(a), the data points of Xq 1 ^ for q = 4 with N = 96, 144 
and 192 seem to fall on the same straight line as the data points for q = 2 and 3 fall on. 
This suggests that the power law behavior like eq. (16) is obeyed by not only for q = 2 and 3 
but for all values of q > 1 at sufficiently small values of q/N. Moreover, the parameters A^\ 
x^\ Arfl and x^\, which satisfy x£\ < x^ and A^\ < \ l ea d to the intersection of the 
two straight lines corresponding to the power law dependences for q = 1 and q = 2 and 3 at 
q/N ~ 5.7 x 10~ 4 . This suggests that the two apparent power law behaviors for q = 1 and 
q > 1 merge into a single power law behavior at small values of q/N. 

As in the case of the trivial knot, the separation of the power law dependences of for 
q = 1 and q = 2 and 3 appears for the trefoil knot and the separation is attributed to the 
topological constraints. The pattern of the separation for the trefoil knot is, however, different 
from that for the trivial knot. In the case of the trefoil knot, A^ > A^\ and x^ > x£\ 
hold, while A^ < A^\ and x^ < x>\ hold for the trivial knot: In Fig. 2(a), the straight 
line corresponding to the power law dependence for q = 1 lies above that for q = 2 and 3, 
while the straight line for q = 1 is located below that for q > 1 in Fig. 1(a). The separation of 
the power law dependences for the trefoil knot leads to the following characteristic behavior 
of the relaxation rates. As A?" increases, the wave number q — <?min 

which gives the slowest 

relaxation rate Xq ^ for each N changes from q m ' m = 2 for A^ < 96 to q m i n = 1 for N > 144. 
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Because X q given by eq. (12) shows the same crossover in Fig. 2(b), the crossover in X q 
can be considered to correspond to the the crossover in {R q ■ R- q ): The largest (R q ■ R- q ) 
for each N is given by q = 2 for small values of N and q = 1 for large values of N. Because 
{Rq ■ R- q ) roughly corresponds to the mean square of the end-to-end distance of a partial 
chain with N/ {2q) segments as mentioned before, the crossover indicates that the end-to-end 
distance of a partial chain with N/2 segments is shorter than that with N/4 segments for 
small values of N and is longer for large values of N. This change can be seen in Fig. 3, 
which shows the longest end-to-end vector of a partial chain with N/2 segments and that 
with Nf 4 segments in snapshots of the configurations of a ring polymer with the trefoil knot 
for N = 24 and N = 192. The longest end-to-end vector of a partial chain with N/2 segments 
is shorter than that with N/4 segments for N = 24 and is longer for N = 192. In Fig. 3, the 



Fig. 3. Snapshots of the configurations of the ring polymer with the trefoil knot for (a) N — 24 and 
(b) iV = 192. The longest end-to-end vector of a partial chain with N/2 segments and that with 
N/4 segments are shown as solid and broken arrows, respectively. In each figure, the plane of 
perspective projection is parallel to the first and second principal axes of the moment of inertia 
tensor for the four end segments of the two end-to-end vectors. 




(a) 
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ring polymer with the trefoil knot is in a "uniform" state, where the knotted part is expanded 
widely along the polymer, for N = 24 and is in a "phase segregated" state, where the knotted 
part is localized to a part of the ring polymer and the rest of the ring polymer behaves like 
a ring polymer with the trivial knot, for N = 192. 8 ' 14 ' 15 Thus, the crossover from the state 
with q m i n = 2 to that with q m in = 1 corresponds to the localization of the knotted part to a 
part of the ring polymer. The localization of the knotted part is considered to be driven by 
the entropy gained by the unknotted part which surpasses the entropy lost by the knotted 
part. 

The localization of the knotted part of the ring polymer with the trefoil knot suggests that 
the proportion of the knotted part to the remainder part, which behaves like a ring polymer 
with the trivial knot, decreases as ./V increases. Therefore, the behavior of the relaxation rates 
of the ring polymer with the trefoil knot is expected to approach that with the trivial knot as 
N increases. From Figs. 1(a) and 2(a), it can be seen that for the trefoil knot approaches 
that for the trivial knot from above as N increases. In contrast, A^ with q = 2 and 3 for the 
trefoil knot approach X^ with q > 1 for the trivial knot from below. As discussed above, for 
each of the trivial and the trefoil knots, the (//iV-dependence of X^ is expected to show a single 
power law behavior at small values of q/N. Therefore, it is expected that the (//iV-dependence 
of X^p shows the same single power law behavior at small values of q/N independently of the 
knot type for sufficiently large values of N. 

4. Summary and Discussion 

In this paper, the relaxation rates of a single ring polymer with the trivial knot and the 
trefoil knot are studied through Brownian dynamics simulations for various values of the num- 
ber of the segments of the ring polymer N. In the case of a single ring polymer, each relaxation 
mode is associated with a wave number q because the ring polymer has the translational in- 
variance along the polymer chain. The slowest relaxation rate A^ for each wave number q 
is estimated by the least square fit of C q (t) given by eq. (7) to the double exponential decay 
(11). The linearization approximation A ? LA ^ to X^ is also estimated from the static correla- 
tion C q (0) on the basis of eq. (12). The behavior of the distribution of A 9 LA - ) is qualitatively 
similar to that of A^ , although each value of A ? LA ^ is larger than the corresponding value of 
Aq 1 ^. Thus, in the case of a single ring polymer, the linearization approximation is useful for 
studying the behavior of the relaxation rate distribution qualitatively. 

For the trivial knot, the relaxation rate A^ appears to behave as A^\l/N) x< i ) for q = 1 
and A^\(q/N) x >i for q > 1 at small values of q/N, where < A { >\ and < x { >\. The 
exponents and x>\ are similar to that found for a single linear polymer chain. Even in 
the case of the trivial knot, the effect of the topological constraints appears as the separation 
of the power law dependences of A^ on q/N for q = 1 and q > 1. 

In the case of the trefoil knot, the topological constraints causes the separation of the 
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power law dependences of X^ for q = 1 and q = 2 and 3 as in the case of the trivial knot. In 
this case, the amplitudes and the exponents of the power law behaviors A^ ~ A^\l/N) x i ' 
for q = 1 and A^ ~ A^^q/N^ 2 ' 3 for q = 2 and 3 satisfy the relations A^ > A^\ and 
> x^\, which is different from the case of the trivial knot. As the consequence of the 
separation of the power law behaviors, the wave number q = q m i n , which gives the slowest 
relaxation rate X^ for each TV, is given by q m i n = 2 for small values of N and q m i n = 1 
for large values of N. This crossover corresponds to the change of the structure of the ring 
polymer from a state for small N, where the knotted part is extended along the ring polymer, 
to another state for large N, where the knotted part is localized to a part of the ring polymer 
and the rest of the ring polymer behaves like the ring polymer with the trivial knot. It is 
expected from the localization of the knotted part and the estimated parameters of the power 
law behaviors that the separated power law behaviors of the relaxation rates for the trivial 
knot and the trefoil knots eventually merge into a single power law behavior in the limit of 
N — > oo, where the effects of the topological constraints vanish. 

In this paper, we only consider the single ring polymers with the trivial knot and the trefoil 
knot, which are the simplest knots. Even for these simple knots, the effects of the topological 
constraints on the relaxation rate distribution are clearly observed. Therefore, it is interesting 
to study the relaxation rates distribution of single ring polymers with more complicated knots 
or ring polymers with links. The study in this direction is in progress. The topological effects 
seem to vanish in the limit that the number of the segments of the ring polymer goes to 
infinity because of the localization of the knotted part. The study of the localization of the 
knotted part of a single ring polymer is also in progress by using the average structure of the 
ring polymer defined self-consistently as the average structure of all the sampled structures, 
which are translated and rotated to minimize the mean square displacement from the average 
structure itself. 

Acknowledgments 

The authors are grateful to Professor T. Deguchi and Dr. K. Hagita for fruitful discus- 
sions. This work was partially supported by the 21st Century COE Program; Integrative 
Mathematical Sciences: Progress in Mathematics Motivated by Social and Natural Sciences. 



12/13 



J. Phys. Soc. Jpn. Full Paper 

References 

1) P. G. de Gennes: Scaling Concepts in Polymer Physics (Cornell University Press, Ithaca, 1984). 

2) M. Doi and S. F. Edwards: The Theory of Polymers Dynamics (Oxford University Press, Oxford, 
1986). 

3) K. Hagita and H. Takano: J. Phys. Soc. Jpn. 71 (2002) 673. 

1) K. Hagita and H. Takano: J. Phys. Soc. Jpn. 72 (2003) 1824 and references therein. 
3) L. H. Kauffman: Knots and Physics (World Scientific, Singapore, 1993). 
3) J. des Cloizeaux: J. Phys. (Paris) Lett. 42 (1981) L433. 
1) J. M. Deutsch: Phys. Rev. E 59 (1999) R2539. 
3) A. Y. Grosberg: Phys. Rev. Lett. 30 (2000) 3858. 
J) T. Dcguchi and K. Tsurusaki: Phys. Rev. E 55 (1997) 6245. 
)) M. K. Shimamura and T. Dcguchi: Phys. Rev. E 64 (2001) 020801(R). 
L) M. K. Shimamura and T. Dcguchi: J. Phys. A 35 (2002) L241. 
1) M. K. Shimamura and T. Deguchi: Phys. Rev. E 65 (2002) 051802. 

3) H. Matsuda, A. Yao, H. Tsukahara, T. Dcguchi, K. Furuta and T. Imai: Phys. Rev. E 68 (2003) 
011102. 

1) E. Orlandini, M. C. Tesi, E. J. Janse van Rensburg and S. G. Whittington: J. Phys. A 31 (1988) 
5953. 

3) E. J. Janse van Rensburg and S. G. Whittington: J. Phys. A 24 (1991) 3935. 
3) S. R. Quake: Phys. Rev. Lett. 73 (1994) 3317. 
1) P.-Y. Lai: Phys. Rev. E 66 (2002) 021805. 

3) P.-Y. Lai, Y.-J. Sheng and H.-K. Tsao: Phys. Rev. Lett. 87 (2001) 175503. 
}) S. A. Wasserman and N. R. Cozzarelli: Science 232 (1986) 951. 
)) S. Y. Shaw and J. C. Wang: Science 260 (1993) 533. 

L) A. Stasiak, V. Katritch, J. Bcdnar, D. Michoud and J. Dubochct: Nature 384 (1996) 122. 
I) H. Takano and S. Miyashita: J. Phys. Soc. Jpn. 64 (1995) 3688. 
3) S. Koseki, H. Hirao and H. Takano: J. Phys. Soc. Jpn. 66 (1997) 1631. 
1) H. Hirao, S. Koseki and H. Takano: J. Phys. Soc. Jpn. 66 (1997) 3399. 
3) K. Hagita and H. Takano: J. Phys. Soc. Jpn. 68 (1999) 401. 
3) K. Hagita, S. Koseki and H. Takano: J. Phys. Soc. Jpn. 68 (1999) 2144. 

1) K. Hagita, D. Ishizuka and H. Takano: J. Phys. Soc. Jpn. 70 (2001) 2897. 
3) P. E. Rouse, Jr.: J. Chcm. Phys. 21 (1953) 1272. 
)) B. H. Zimm: J. Chcm. Phys. 24 (1956) 269. 
)) G. S. Grest and K. Krcmer: Phys. Rev. A 33 (1986) 3628. 
L) K. Krcmer and G. S. Grest: J. Chcm. Phys. 92 (1990) 5057. 

2) K. Binder: Monte Carlo and Molecular Dynamics Simulations in Polymer Science (Oxford Uni- 
versity Press, Oxford, 1995). 

33) V. Bloomfield and B. H. Zimm: J. Chcm. Phys. 44 (1966) 315. 



13/13 



